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ABSTRACT 

We report the discovery of two extremely low luminosity globular clusters in the Milky Way Halo. 
These objects were detected in the Sloan Digital Sky Survey Data Release 5 and confirmed with 
deeper imaging at the Calar Alto Observatory. The clusters, Koposov 1 and Koposov 2, are located 
at ~ 40 — 50 kpc and appear to have old stellar populations and luminosities of only My ~ —lmag. 
Their observed sizes of ~ 3 pc are well within the expected tidal limit of ~10 pc at that distance. 
Together with Palomar 1 , AM 4 and Whiting 1 , these new clusters are the lowest luminosity globulars 
orbiting the Milky Way, with Koposov 2 the most extreme. Koposov 1 appears to lie close to distant 
branch of the Sagittarius stream. The half-mass relaxation times of Koposov 1 and 2 are only ~ 70 
and ~ 55 Myr respectively (2 orders of magnitude shorter than the age of the stellar populations), 
so it would seem that they have undergone drastic mass segregation. Since they do not appear to 
be very concentrated, their evaporation timescales may be as low as ~ O.ltnubbie- These discoveries 
show that the structural parameter space of globular clusters in the Milky Way halo is not yet fully 
explored. They also add, through their short remaining survival times, significant direct evidence for 
a once much larger population of globular clusters. 
Subject headings: Galaxy:halo - Galaxy: globular clusters 



1. INTRODUCTION 

The population of globular clusters around the Milky 
Way has been studied extensively and the current census 
finds the majority at low latitudes in the inner Galaxy 
(Rgc < 20 kpc). Globular clusters are almost univer- 
sally "old" [i ago « 0.5 1 x t Hubble]; show no con- 
vincing evidence for dark matter, and have character- 
istic luminosities of 10 5 L Q My ~ —8) and typical sizes 
of 3 pc. Yet, the observed range of structural proper- 
ties (e.g. mass, size, and concentration) is quite wide. 
This range is of great interest, as it appears to be de- 
termined by a set of astrophysical processes: the ini- 
tial structure and orbit; subsequent external processes, 
such as galactic tides and dynamical friction; and ensu- 
ing mass segregation, evaporat i on, and core collapse (see , 
e.g., iGnedin fc Ostrikeri Il997l : iMevlan fe Heggiel I1997D . 
Indeed, there has long been a sense that the observed 
population of Galactic globular clusters mainly reflects 
the subset of objects that could survive for ~ tnubbic- In 
individual cases, there is clear evidence for intern al re- 
shaping processes (as in M15 , ISosin fc K ing 199 7]) an d 
tidal disruption (as in Pal5, Odcnkirch en et aTl 12001). 
Within this context, identification and study of globular 
clusters with extreme properties is undoubtely of great 
interest. 

Our census of objects at the outskirts of the Milky 
Way has increased rapidly in the last few years, mostly 
based on large-area CCD surveys such as the Sloan 
Digital Sky Survey (SDSS; lYork et all I2000L Recent 
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added two faint and extended objects that may be Milky 
Way globular clusters. The newcomers, Willman 1 and 
Segue 1, both have distorted irregular isopleths, perhaps 
indicating ongoing tidal disruption. Willman 1 seems 
to show so me evidence for da rk matter and metallic- 
ity spread (]Martin et al.l l2007h . casting some doubt on 
whether it is a globular cluster at all. 

Here, we announce the discovery of two new, distant, 
extremely faint and compact (~ 3 pc) globular clusters, 
named Koposov 1 and Koposov 2, first detected in SDSS 
Data Release 5 (DR5) and subsequently confirmed with 
deeper imaging at Calar Alto. The total luminosity of 
Koposov 2 appears to be ~ —1 mag, lower than that 
of the faint est Galactic globular kn own to date, AM 4 
(—1.4 mag. Unman fe Carnevl fT987h . Koposov 1 is not 
much brighter; at My ttot ~ —2 mag, it has the third- 
lowest luminosity. In total, only 3 out of the previously 
known ~ 160 Galactic clusters, have comparably low lu- 
m inosity and small siz es: AM 4, Palomar 1, and Whiting 
1 (| Whiting et alJl2002f ). Willman 1 and Segue 1 also have 
extremely low luminosities, but are an order of magni- 
tude larger. 

In this paper, we describe the deep follow-up data con- 
firming the discoveries and give estimates of the struc- 
tural parameters of the new objects (see Table [TJ. We 
argue that the discovery of these two low-mass globulars 
in less than 1 /5 of the sky may mean that a substantial 
population of such clusters lurks in the outer halo of the 
Milky Way. 

2. DISCOVERY AND OBSERVATIONS 

The two new globular clusters were originally selected 
among other candidates in the course of our systematic 
search for small-scale substructure in the Milky Way 
halo. The aim of the search was to detect all signif- 
icant small-scale stellar overdensities above the slowly 
varying Galactic background that are likely to be ei- 
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Fig. 1. — 3' X 3' SDSS cutout images of Koposov 1 and 2. The bright star in the center of Kop osov 1 is a foregroun d star with V ~ 14. 5 m 
and large proper motion (/tcoAM) ~ ( — 32, — 12) fias yr~ 1 , according to the USNO-B1 catalog (Monet ct al. 2003). The bright extended 
object near the center of Koposov 2 is a background galaxy. 



TABLE 1 

Properties of Koposov 1 and Koposov 2 



Parameter 



Koposov 1 



Koposov 2 



Coordinates (J2000) 11:59:18.5 +12:15:36 07:58:17.00 +26:15:18 

Coordinates (V, 6) (260.98° , 70.75°) (195.11° , 25.55°) 

Distance ~ 50 kpc ~ 40 kpc 

Size ~ 3 pc ~ 3 pc 

M v ~ -2 m l m 

Relaxation Time ~ 70 Myrs ~ 55 Myrs 

Tidal radius ~ 11 pc ~ 9 pc 



ther dwarf spheroidal galaxies or globular clusters. A 
detailed description of the algorithm and its efficiency 
will be provided in a future paper, and we only present 
here a brief outline of the method. The algorithm is 
based on the so-called Difference of Gaussians method , 
first developed in Computer Vision (|Babaud et al.lH98ol 
iLindenberd Il998h . Starting from a flux-limited cata- 
log of stellar positions, the number-counts map in (a, 8) 
plane is convolved with a filter optimized for the detec- 
tion of overdensities, na mely the difference o f two two- 
dimensional Gaussians (|Koposov et al. 2007). Having 
zero integral, the kernel guarantees that the convolution 
with a constant (or slowly varying) background will re- 
sult in zero signal. When the data contain an overdensity 
with a size comparable to the size of the inner Gaussian, 
the filter will be close to optimal. 

We applied this filtering procedure to the entire stellar 
subset of the DR5 source catalog with r < 22 m , g — r < 
1.2 m . In our analysis we used the photometry cleaned 
by switching on quality flags as described in SDSS SQL 
pages 3 This minimizes the influence of various artefacts 
including those caused by proximity of very bright or 
extended objects. In the resulting map that had been 
convolved with a 2' kernel, we found two very compact 
objects among other overdensities ranked highly accord- 

l! http: //cas . sdss . org/dr5/en/help/docs/real query . asp#f lags 



ing to their statistical significance. Figure [T] shows the 
SDSS images, and Figure [2] shows the spatial distribu- 
tion of extracted sources, where central concentrations 
of stars are clearly visible. These concentrations are de- 
tected at high level of significance. The areas of 1' radius 
marked by circles centered on Koposov 1 and 2 plotted in 
Figure [5] contain 22 objects and 23 objects, respectively, 
while mean density of g-r< 0.6, r> 20mag stars should 
produce approximately 2.5 objects, which implies a high 
statistical significance of the overdensities; for pure Pois- 
son distribution of objects, the probability to find such 
group of stars in all DR5 is around 10~ 9 . 

The differential Hess diagrams for stars within 2'. 5 ra- 
dius centered on the objects are shown in Figured There 
is a clear excess of blue stars (g — r < 0.5), which we in- 
terpret as main-sequence turnoff stars at r ~ 22, which 
roughly corresponds to distances of ~ 50kpc. 

To confirm the nature of discovered candidates and 
quantify their structural and population properties, we 
acquired follow-up GTO observations in 2007 January 
on the 2.2m telescope at Calar Alto using the CAFOS 
camera. This camera has a 2k x 2k CCD with a 16' x 16' 
field of view and a pixel scale of 0" .5pixel~l. We ob- 
served each object for a total of 2 hr in Johnson B 
and 1.5 hr in Cousins R. The integrations were split 
into five individual dithered exposures for cosmic ray 
and bad pixel rejection. The observations were car- 
ried out in good photometric conditions with a seeing of 
1" — 1.3". The data were bias-subtracted and flat-fielded. 
The individual frames were WCS-aligned, drizzled, and 
median-combined using our softwar e and the SCAMP 
and SWARP programs (|Bertml 12006). The combined B 
band images of the objects are shown in Figure 31 

The central stellar overdensities are clearly corrobo- 
rated by the Calar Alto photometry, which is nearly 2 
mag deeper than the original SDSS data. While the 
follow-up data are quite deep, the stars are subject to sig- 
nificant crowding, due to the compactness of the clusters. 
Therefore, for the purposes of robust source detection 
and photometry, we used the DAOPHOT/ALLSTAR 
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Fig. 2. — The spatial distribution of the objects in the area of Koposov 1 and Koposov 2. All objects classified as stars with colors 
g-r)< 0.6 m and r> 20 m in the area 0.3° X 0.3° are shown. The circles with 1' radii centered on the objects are overplotted. 
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Fig. 3. — The residual g — r versus g Hess diagrams of the clusters from the SDSS data. In each case the residual Hess diagram is 
costructed by subtracting the normalized background Hess diagram from the Hess diagram of stars lying within 2'. 5 radius from the centers 
of objects 



software (|Stetsonl fl987) . To get the absolute calibra- 
tion of the photometry from each frame, we cross- 
matched the DAOPHOT sources with the SDSS cata- 
log using the Virtual Obser vatory resource SAI CAS 4 
(|Koposov & Bartunovf2 006). To convert the Sloan g and 
r magnitudes into the Johnson-Cousins pho tometric sys- 
tem, w e used the conversion coefficients from lSmith et all 
( 2002). The resulting B versus B — R color-magnitude 
diagrams (CMDs) of the central regions of the objects 
together with the CMDs of the comparison fields, ex- 
tending to B ~ 23. 5 m — 24 m , are shown in the Figure 
The median photometric accuracy of the data is 0.05-0.1 

4 |htt p : //vo . astron et ■ ru/cas] 



mag. The CMDs clearly show the presence of the main 
sequences near the centers of the objects, while they are 
absent in the the comparison fields. The statistical sig- 
nificance of the overdensities is also clearly supported by 
the new data. The CMD of objects within 2' from the 
center of Koposov 1 contain 96 objects, while the back- 
ground density inferred from the comparison field should 
give around 23 objects, which gives a 15 a deviation. For 
Koposov 2 , the number of objects within 1.2' is 92, while 
the background density from the comparison field should 
produce around 24 objects, which gives a 14 sigma devi- 
ation. In the next section we will discuss the properties 
of the objects which can be derived from the follow-up 
data. 
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Fig. 4. — B band Calar Alto view of Koposov 1 and Koposov 2. The 2' X 2' images are centered on the clusters (north is up, cast is left). 
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Fig. 5. — Left panel, left half : B vs. B — R CMDs derived from the Calar Alto data for stars lying within 2' of Koposov 1 with 8 Gyr and 
[Fe/H]=-2 Girardi ct al. (2000) isochrones overplotted. Left panel, right half. For comparison, the CMDs of stars in the annulus centered 
on Koposov 1 defined by ra dii 3.2' and 3. T are p lotted. Right panel, left half. B vs. B — R CMDs of stars lying within 1.'2 of Koposov 
2 with 8 Gyr and [Fc/H]=-2 Girardi ct al. (2000) isochrones overplotted. Right panel, right half, for comparison, the CMDs of stars in the 
annulus centered on Koposov 2 defined by radii 2' and 2.3'are plotted. 



3. PROPERTIES 

The CMDs of the objects from the Calar Alto data 
(Fig. [5]) clearly show a distribution of stars which can be 
attributed quite convincingly to an old main sequence. 
In the case of Koposov 1, the main-sequence turnoff is 
clear-cut, while for the second cluster it is not so well 
defined. To estimate the distances to the objects, we 
overp lot in Figure [ j] the 8 Gyr [Fe/H]=-2 isochrones 
from I Girardi et all (|2000D . For Koposov 1, this gives 
a distance of 50 kpc. For Koposov 2, the estimate is 
40±5 kpc, but it is not well constrained due to a lack of 
main-sequence turnoff stars. The angular diameters of 
the clusters are < 0.5', which translates into a physical 
size of r ~ 5pc. Unfortunately, the number of stars de- 
tected in the central regions is not enough to precisely 
measure half-light radii of the objects; our best estimate 
is rh ~ 3pc. For Koposov 1, we subtracted the bright 
foreground star near the center, integrated the light of 
the whole cluster in apertures and fitted it to a Plummer 
profile with r/j = 3pc. For Koposov 2, we performed a 
maximum likelihood fit with ~ 3pc. Moreover, the 
minuscule number of stars in both clusters does not al- 
low us to firmly establish their total luminosities. Our 



estimate of — 1 > My > —2 is based on the absence 
of the giants in these clusters and the visible similarity 
of the CMDs to that of t he lowest luminosity glo bular 
cluster AM4 {M v = -1.6. Unman fc Carnevill987| ). We 
checked that estimate by a simple Monte Carlo exper- 
iment: using the Salpeter IMF and Girardi isochrones 
we simulated fake clusters and deduced that the clusters 
with —1 > My > —2 have a number of stars within 1.5-2 
mag below the turnoff is close to the observed number of 
stars (50-70) in our objects. We must say also that due 
to the intrinsic faintness of the clusters and low num- 
ber of stars in them the estimates of the total luminosity 
and especially the age have large uncertainties. However, 
with the existing data we cannot do any better. Much 
deeper and more accurate photometry may be required 
to get precise age/luminosity measures. The spectro- 
scopic observations would be interesting in constraining 
the metallicity of these objects, which is currently com- 
pletely unknown. 

We note that the CMD of Koposov 1 shows several 
stars brighter and bluer than the tentative main-sequence 
turnoff, which we interpret as blue stragglers. This hy- 
pothesis is not implausible considering the low luminos- 
ity of the cluster and taking into account the observed 
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Fig. 6. — Right ascension v s. distance for the A a nd C branches 
of the Sagittarius stream (sec Bclokurov ct al. 2006j). The position 
of Koposov 1 is marked by a star. 

anti-correlation between the frequency of bl ue stragglers 
and t he luminosity of the globular cluster (|Piotto et al.l 
l200l . 

The distance and the position of Koposov 1 suggest 
that this cluster may be related to the Sagittarius tidal 
stream. Its location is a good match to t he distant 
tidal arm discovered in lBelokurov et all ((2006). Figure [6] 
shows the arms of the Sagittarius stream in the DR5 slice 
around S ~ 10° and the position of Koposov 1. 

4. DISCUSSION 

Figure [7] shows Koposov 1 and 2 on the size- luminosity 
plane along with other Galactic globular clusters. This 
illustrates how unusual Koposov 1 and 2 are in their 
structural properties. It appears that the detection of 
these clusters contributes to growing evidence for a large 
population of small and extremely faint objects (includ- 
ing Palomar 1, AM 4, E3 and Whiting 1). There is a 
clear indication as well that this sub-population of glob- 
ular clusters may have significantly yo unger ages than 
classi cal globular clus ters: Palomar 1 (ISaraiedini et al.l 
120061 ) and Whiting 1 (jCarraro et aT1l2007l ) have ages be- 
tween 4 and 6 Gyrs. The current estimate for the age 
of Koposov 1 is « 8 Gyr, and the age of E3 globular 
cluster is « 10 Gyr. This group of clusters is also quite 
apparent on the galactocentric distance versus luminos- 
ity plane shown in Figure [H At least 2 out of these 5 
unusual clusters (Whiting 1 and Koposov 1) seem to be 
associated with the Sagittarius dwarf galaxy. Two quan- 
tities that are crucial for the long term evolution and 
survival of Koposov 1 and 2 are the relaxation time and 
the expected tidal radius. For the ha lf-mass r elaxat ion 
time, we find usin g equation (2-63 ) of ISpitzerl (|1987f ) or 
equation (72) of iMevlan k. Heggid (|1997h . 
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respectively for Koposov 1 and Koposov 2. Here, we 
have assumed L « 200L Q , M/L « 1.5, (m*) ~ 0.6M Q 
and N=500 for Koposov 2, while for Koposov 1, we have 
assumed twice as many stars, using the observational 
estimates of § [3] This means that both clusters have ex- 
tremely short relaxation times, less than 1% of imibbie 
and t r h w 0.01t age ». The most immediate effect of two- 
body relaxation is mass segregation, which should be 
quite drastic given the apparent stellar population age. 



Fig. 7. — Size vs. absolut e magn i tude p lot for Galactic globular 
clusters. The data from the Harris ( 1996) catalog are plotted with 
diamonds. Squares mark the locations of the recently discovered 
globular clusters Willman 1, Segue 1 and Whiting 1. Koposov 1 
and 2 are shown as stars. 

The expecte d tidal truncation o f these clusters occurs 
at (see, e.g.. llnnanen et afll!983[ ) 
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where we have assumed an orbital eccentricity of 0.5, 
and that the clusters are now near apocenter (hence 
Rperi ~ 16 kpc), a Milky Way circular speed of 190 
kms^ 1 at 16 kpc and a cluster (stellar) mass of 600 and 
300 Mq for Koposov 1 and 2, respectively. Hence, the 
detectable extent of the globular clusters (3 pc) falls well 
within the tidal limit. From this argument, the clusters 
are under no threat of destruction by tidal forces. Al- 
though formal profile fits are not feasible with so few 
stars, the stellar distributions (see Figs. Q] and [3| are 
well localized, but not centrally concentrated by globu- 
lar cluster standards; a core to tidal radius ratio of the 
observed stellar distribution of 4 seems reasonable, im- 
plying a concentration parameter of c = log(r t /r c ) w 0.5. 
For such low concentrations, the evaporation timescale 
t cv , which is the time-scale over which two-body re- 
laxation drives stars to beyond the escape velocity, is 
t cv « 1.5t cc w 12t r h (where t cc is core collapse time) (Fig- 
ure 17 and 19 in iGnedin et al.lfl999T )For Koposov 1 and 
2, this implies evaporation time-scales of 0.7 and 1.1 Gyr, 
respectively. This estimate of t ev ~ O.lf. Hubble may be an 
underestimate, if the brightest stars which we observe 
are more concentrated than the faint stars due to mass 
segregation; then the total mass and half-mass radius 
can be larger. Nonetheless, this estimate makes it clear 
that the present structural and dynamical state cannot 
have prevailed, even approximately, for a time-span of 
~ 10 Gyr. The above arguments hold irrespective of 
whether Koposov 1 and Koposov 2 were once part of a 
satellite galaxy, because they are mostly derived from in- 
ternal evolution factors. This discrepancy of time-scales 
is more pronounced in Koposov 1 and 2, because their 
relaxation time-scales are shorter than those of Palomar 
1 and Whiting 1, which in any case have accurate pho- 
tometry suggesting younger ages of ~ 4 — 6 Gyr. 

At face value, Koposov 1 and 2 have survival times in 
their current state of ~ 0.1<Hubbie, and were found in a 
search of 20% of the whole sky (SDSS DR5). The naive 
multiplication of these factors points to a large parent 
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Fig. 8. — Galactocentric distance vs. magnitude plot for Galactic 
globular clusters. Symbols are as in Fig. [7] 

population of ~ 100 objects. The most likely reservoir 
for this parent population is the globular clusters, and 
possibly even old open clusters, in satellite galaxies that 
have been accreted, like the Sagittarius. In objects like 
Koposov 1 and 2, it is clear that the very short relaxation 
and evaporation times must lead to drastic mass segrega- 
tion and the expulsion of basically all low-mass stars (this 
line of reasoning lead us to the modest M/L w 1.5). This 
gives new life to the view that truly many of the accreted 
globular clusters must have been destroyed. Yet, it is also 
clear that the actual dynamical prehistory and future of 
these clusters requires much more careful modelling. The 
small number of stars makes them ideal subjects of direct 
iV-body calculations. But regardless of their dynamical 
evolution, these clusters manifestly demonstrate that the 
parameter space of globular clusters in the Milky Way is 
not yet fully explored. 
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